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ABSTRACT 

In this letter we present, for the first time, a study of star formation rate, gas fraction and galaxy 
morphology of a constrained simulation of the Milky Way (MW) and Andromeda (M31) galaxies, 
compared to other MW-mass galaxies. By combining with unconstrained simulations we cover a 
sufficient volume to compare these galaxies environmental densities ranging from the field to that of 
the Local Group (LG). This is particularly relevant as it has been shown that, quite generally, galaxy 
properties depend intimately upon their environment, most prominently when galaxies in clusters are 
compared to those in the field. For galaxies in loose groups such as the LG, however, environmental 
effects have been less clear. We consider the galaxy’s environmental density in spheres of 1200 kpc 
(comoving) and find that whilst environment does not appear to directly affect morphology, there is 
a positive trend with star formation rates. This enhancement in star formation occurs systematically 
for galaxies in higher density environments, regardless whether they are part of the LG or in filaments. 

Our simulations suggest that the richer environment at Mpc-scales may help replenish the star-forming 
gas, allowing higher specific star formation rates in galaxies such as the MW. 

Subject headings: galaxies: formation - galaxies: evolution - cosmology: theory - methods: numerical 


1. INTRODUCTION 

During the last decades, observational studies found 
indications that the environment where galaxies form 
plays an important role on the determination of their fi¬ 
nal properties. For example, elliptic al galaxies are more 


strongly clustered than spirals (e.g . Dressier 1980 Her- 
mit et ak1|1996 Guzzo et al. 1997); and correlations be- 


tween environment and the photometric colour or lumi¬ 


nosity of a galaxy have also been measured (Blanton et al. 
20031. The suppression of star formation in clustered 


enviro nments also extends to groups (Lewis et al. 2002 
Coenda et al.|2014 ) and the morphology trends, i.e. that 
early type galaxies lie closer t o halo centres, exten d as 
a continuum to group scales (Girardi et al. 20031, al¬ 
though at stellar mas ses closer to that o f the MW this 
is more controversial (Ziparo et al. 2013) and some au¬ 
thors suggest large scale features such as fil aments may 


be more important than the group (e.g. Balre et al. 


20131. Group galaxies also app ear deficient i n cold gas 


] Bess & Wilcots 2013), although McGee et al. (2008) find 
that whilst the disc fraction depends on the group envi¬ 
ronment, the other properties of the discs appear largely 
unaffected. 

On the other hand, theoretical studies in the context of 
A-Cold Dark Matter (ACDM) show that the merger rate 
of galaxies depends on environment, such that galaxies 


est fractions of mergers (e.g. Maulbetsch et al. 2007 

de 

Ravel et al. 2009; Darg et al. 2010; Ellison et al. |20 

L0; 

Fakhouri el^jiLJ 

2010( |Lin et al. 2010; |Tonnesen & Cen 

2012 Jian et a 

. 2012). Mergers are indeed known to 
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Toomre & Toomre 1972), and could be responsible for 
at least some of the observed trends of galaxy properties 
with environment. 

Our MW lives in a rich group environment, known as 
the “Local Group” (LG). Andromeda (M31), the other 
large spiral in the LG, has a similar mass to the MW, 
and lies at less than a Mpc away. Such spiral galaxies are 
prevalent in the Universe, and thus a fundamental test 
of the ACDM cosmol ogical fr a mework is that neither the 
local universe (e.g. |Nuza et al. |2014a| nor the MW and 
M31 should be extremely improbable objects. 

The numerical simulation of MW-like galaxies is an 
active ar ea of research, with zoom-in techniques (Katz 
& White 1993) allowing the collisionless and collisional 


components of cosmological density fluctuations to be 
followed over almost the full 13.7 billion years of cos¬ 
mic time. Modern computers enable simulations that 
describe the internal properties of galaxies so that time 
variations due to mergers, interactions and infall can in¬ 
deed be quantified using an acceptable number of reso¬ 
lution elements. 

The precise initial conditions (ICs) to produce late type 
galaxies such as the MW and M31 are not believed to be 
strict as spiral galaxies are the most abundant type in 
our Universe. In most modern studies (e.g. Scanna¬ 


pieco et aL|2009HGuedes et al.||201l]|Aumer et aE 2013 


Stinson et al.|2013 Vogelsberger et ai.||2014[ ) the ICs are 

chosen as halos of MW mass and isolated from other mas¬ 
sive halos that can destroy discs at late times. Indeed the 
primary uncertainty of such simulations lies in the unre¬ 
solved baryonic phys ics that res ult in star formation and 
feedback (see Creasey et al. 2011), and different codes ap¬ 
plied to the same halo ca n produce substantial var iation 
in the stellar component (Scannapieco et al. 2012). 

Whilst the interaction of the 1V131-MW pair is not yet 
expected to affect internal galaxy properties as they are 
still too distant, their environments on Mpc scales have 






















































































2 


merged, and this overdense environment may affect the 
pair. One method to test such effects is to identify anal¬ 
ogous pairs produced by cosmological ICs, and com pare 
to mo re isol ated systems. Following this approach, |Few| 
et al. (20121 found that for simulated MW-mass galax- 
ies the difference between the held and loose groups ap¬ 
pears marginal, and in fact they detected no visible differ¬ 
ence between MW-likc held ga laxies a nd th ose that reside 
in loose groups. Garrison-Kimmel et al. (2014) find no 
difference in concentrations or stellar masses within the 
virial radius, but do hnd an in creased num ber of ‘back- 
splash’ galaxies (see e.g. Knebe et al.||2011) - those that 
have escaped the virial radius - lor pairs such as the MW 
and M31. 

Although those simulations allow assessment of the ef¬ 
fects of environment on the formation of MW-mass galax¬ 
ies, they do not fully exploit the detailed constraints we 
have from the present-day dynamics of our local uni¬ 
verse. One method to ameliorate this is to utilise the 
power of constrained ICs that reproduce the observed 
dynamical properties of the LG in combination with hy- 
drodynamical simulations, which has not previously been 
attempted. This provides a pertinent test to quantify the 
influence of environment, as well as allowing us to dis¬ 
cern directly the physical processes in play during the 
formation of galaxies in a LG-like environment. 

In this letter, we present results from a simulation con¬ 
strained both to form a LG analogue and to match the 
velocity field of the local universe. This allows us to con¬ 
trast the properties of the simulated LG galaxies with 
other galaxies of similar stellar mass. For our LG sim¬ 
ulation we have taken (a version with slightly greater 
co verage o f) t he ICs used to produce a MW-M31 pair 


Nuza et al. (2014b) and Scannapieco et al. (submit¬ 


ted) in a volume resembling the distribution of matter 
of our local Universe. To enlar ge and contrast ou r sam- 


ple we include the 8 galaxies of Scannapieco et al. (2009 
S09 hereafter) which use the same code. The latter sim¬ 
ulations are themselves resim ul ations o f the 8 halos of 
the Aquarius Project ((Springel et al.p008j) that result in 
isolated galaxies at z = 0. 

This letter is organised as follows. In Section [2] we 
describe the simulation code that was used to evolve the 
ICs of both the LG and Aquarius, and the sample of MW- 
analogues that are produced. In Section [3] we analyse 
the evolution of the environment of these galaxies and 
attempt to discern its effect on basic galaxy quantities 
such as stellar mass, gas mass and star formation rate. 
In Section 0 we discuss and conclude. 

2. SIMULATIONS 

2.1. Initial conditions 

The ICs used for our LG simulations are part of 
the CLUES (Constrained Local UniversE Simulations 3 ) 
Project. The ICs reproduce, by construction, the known 
dynamical properties of our local environ ment at z = 0 
(jGottlober et al.||2010| |Yepes et al.||2014|) , and are con¬ 
sistent with a ACDM universe with W MAP-5 parame¬ 
ters: Um = 0.279 (matter density), Ua = 0.721 (dark 
energy density), flbar = 0.046 (baryon density), H 0 = 
100/ikms -1 Mpc -1 with h = 0.7 (Hubble parameter), 
and <7g = 0.8 (normalization of the power spectrum). 


http://www.clues-proj ect.org/ 


The ICs use the zoom-in technique, where particles 
are placed in a periodic box of 91 Mpc on a side and our 
high-resolution region is of 10 Mpc radius at z = 0, with 
a mass resolution of 2.8x 10 6 M 0 and 5.6x 10 5 M 0 in dark 
matter and gas particles respectively, and a gravitational 
softening length of 0.7 kpc, fixed in physical coordinates 
since z = 3 and fixed in comoving coordinates at earlier 
ti mes. The IC phases are the same as for |Nuza et al. 


(2014b), but the zoom region is slightly larger, causing 


some variation in the simulated evolution. 

The Aquarius simulations that we use to compare with 
our LG galaxies are the eight galaxies first presented 
in S09. These are the hydrodynamical counterparts of 
the galaxies of the Aquarius Project, selected to have 
formed in isolated environments by requiring that they 
have no neighbour exceeding half their mass within a 
sphere of 1.4 Mpc at z = 0. The cosmological parame¬ 
ters are slightly different to those used in our LG simu¬ 
lation: U M = 0.25, U A = 0.75, U bar = 0.04, h = 0.73, 
and (t 8 = 0.9. The mass resolution and softening lengths 
adopted are however similar in the two samples, and we 
have used the same set of input parameters for star for¬ 
mation and feedback. The different cosmological param¬ 
eters used are not expected to compromise our results as 
the adjustment to the cosmology is rather minor in com¬ 
parison with the dispersion in the evolutionary trends. 

2.2. Simulation code 

We use an extended version of the Tree-PM Smoothed 


Particle Hy drodynamics (SPH) code GADGET3 (Springel 
et al.|2008) that includes metal-dependent cooling, chem¬ 


ical enrichment and feedback from Type II and la super- 


ground field (Haardt & Madau 1996 

). The model has 

been developed in Scannapieco et al. 

,2005 2006) and is 

the same code as lor jvluza et al. (j2oi 
In previous work, we have shown th 
to reproduce the formation of galaxy 

4li|). 

it our model is able 
discs from cosmo- 


the angular momentum problem. Some limitations ex¬ 
ist, such as simulated galaxies tend to have overl y mas- 
sive bulges (though see the implementation of Aumer 


|et al.||20l3| where the bulge mass is significantly reduced 
via early stellar feedback), however discs do have real¬ 


istic sizes and angular momen tum content (Scannapieco 
et al. 2008 2009 2010| 2012), allowing studies of their 
formation and evolution m relation to those structures on 
larger scales. In a companion paper (Scannapieco et al., 
submitted) we analyse the evolution of the distribution 
of the stellar component, to which we refer the reader for 
details of morphology classification. 


2.3. The Galaxy Sample 

We choose a sample of galaxies with stellar masses in 
the range 10 10 ' 5-115 M 0 , i.e. similar to the MW. For 
our LG simulation this includes the M31 and MW can¬ 
didates 4 G1 and G2 along with four other galaxies that 
we refer to as G3-G6. These galaxies have no ‘contamina¬ 
tion’ within 1200 comoving kpc (ckpc) from their centres, 
i.e. they are sufficiently embedded in the zoom region to 
avoid low-resolution particles. The environment of these 

4 Note that in|Nuza et al. (2014b[i these galaxies were referred 
to as M31 c and \TU ' respectively. 
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Fig. 1.— Column density of gas (blue) overlaid with star forming 
gas ( orange-white ) for the zoom region at z = 0. Red circles delimit 
the R 200 f° r each galaxy of a similar stellar mass to the MW (Gl- 
G6), and the labelled insets enlarge these regions. 1200 kpc is the 
radius of our environment measure, and the width of the image is 
20 Mpc. 


TABLE 1 

Main properties of the LG and Aquarius galaxies at z = 0: 

VIRIAL RADIUS (R20o)i VIRIAL MASS (M 200 ), AND MASSES IN 
STARS, GAS AND STAR-FORMING GAS (M s tar, M gaH AND Mgp), 
WITHIN THE VIRIAL RADIUS. 



R 200 

M 200 

M s t a r 

M gaa 

Msf 


[kpc] 

[ 10 10 Mq] 

[10 10 M 0 ] 

[ 10 10 M 0 ] 

[10 10 M 0 ] 

Gl 

245 

168 

8.1 

6.6 

0.479 

G2 

219 

120 

6.4 

5.2 

0.369 

G3 

211 

108 

6.8 

3.6 

0.078 

G4 

166 

52 

3.3 

1.8 

0.146 

G5 

177 

63 

4.4 

1.8 

0.151 

G6 

174 

60 

3.7 

1.5 

0.090 

Aq-A 

232 

149 

9.2 

4.9 

0.153 

Aq-B 

181 

71 

4.0 

1.7 

0.033 

Aq-C 

237 

161 

11.0 

3.8 

0.119 

Aq-D 

233 

149 

8.4 

3.5 

0.005 

Aq-E 

206 

108 

8.4 

2.6 

0.048 

Aq-F 

196 

91 

7.7 

1.8 

0.012 

Aq-G 

180 

68 

4.5 

1.6 

0.061 

Aq-H 

182 

74 

6.5 

0.6 

0.011 


galaxies in gas and star-forming gas is projected in Fig. l] 
along with insets to enlarge the properties within the 
virial radius. All the galaxies in the LG sample are well 
separated (> 3Mpc), except the G1-G2 pair at 770 kpc. 5 

The eight Aquarius galaxies (denoted AqA-H) inhabit 
disjoint environments and were resimulated individually. 

5 A visualisation of a revolution of this structure and the envi¬ 
ronment selections is available here 


Images of these galaxies can be seen in S09 and for previ¬ 
ous studies of thei r formation historie s and properties w e 
refer the reader to Scannapieco et al. (2009 2010| 2011). 
Table [l] summarises the main properties or the LG and 
Aquarius galaxies at z = 0 inside their respective virial 
radii R 20 O) where the density is 200 times the critical 
density p c at the corresponding redshift. 

The LG galaxies have stellar components that all show 
some measure of rotational support, however when we 
examine the components of angular momenta (e.g. as 
in S09) we see that only half, G2, G3 and G4, have 
a significant fraction (> 17%) of the stellar mass in 
a rotationally-supported disk. Nevertheless, all these 
galaxies exhibit extended gas discs and star formation 
at redshift zero. In the case of the Aquarius simulations 
each galaxy was able to grow extended discs during their 
evolution, but again only half can survive until z = 0, Aq- 
C, Aq-D, Aq-E and Aq-G. In both cases we find the sur¬ 
vival or destruction of discs in the S09 sample was found 
to depend primarily on the occurrence of major mergers 
and the alignment between the angular momenta of the 
stars and gas in the inner regions (see also Scannapieco 
et al., submitted). 


3. ENVIRONMENTAL EFFECTS 

The constrained nature of our simulation allows us to 
explore the possibility that environment plays a role in 
the determination of the properties of galaxies like our 
MW. In this section, we investigate this by comparing the 
properties of the LG and Aquarius galaxies. We compare 
properties that are expected to depend on environment 
(gas mass, stellar mass, star formation rate), to a mea¬ 
sure of the environment. We note that the number of 
galaxies in our samples is too small to determine quan¬ 
titative trends, but that systematic differences may still 
be visible within the sample. 

For our environmental measure we use the ratio of the 
mean density of matter within 1200 ckpc to the mean 
density of matter of the universe, i.e. 


<5i200 = 


(Pm) 

r<1200 ckpc 

^m(z)Pc(2) 


( 1 ) 


To give some physical reference, this radius corresponds 
to how far a baryon travelling at lOOkms -1 traverses in 
12 Gyr (ignoring cosmological expansion), approximately 
identifying the cosmological neighbourhood which can af¬ 
fect a MW-mass halo. This scale was also used a s proxy 
for the 'local volume’ by Garrison-Kimmel et al. (2014), 
and in the case of the LCsimulation is large enough to 
include the halos of both G1 and G2 at the present time, 
but for all other galaxies forms disjoint volumes. We have 
considered scales from 600-1500 ckpc and found this to 
give the least stochastic results. 

Fig- @ compares the environments of our galaxies as 
a function of redshift, giving us a measure of the en¬ 
vironmental ‘assembly’, although we would stress that 
at these large radii we are not seeing halo mergers but 
rather a growth in the richness of the environment. Ap¬ 
parent is that although all the galaxies inhabit similar 
environments at 2 ~ 2-3, the three galaxies Gl, G2 and 
G4 exhibit stronger evolution after z ss 2 and are the 
galaxies with the highest overdensity di 2 oo at z = 0, ap¬ 
proximately twice as overdense as the Aquarius galaxies 


















4 



Fig. 2.— Environmental overdensity 5i200 vs. redshift for the 
galaxies in the LG and Aquarius samples. 


and the remaining LG galaxies. Additionally these three 
all inhabit gaseous filam ents, the fo rmer tw o reside in 
the filament identified by Nuza et al. (2014b) and for the 
latter can be seen in the animated (online) version of 
Fig -0 Such structures h ave previously bee n implicated 
as affecting SF in groups (B ahe et aL| 2013). Since these 
three galaxies recur as outliers later we denote them the 
‘rich sample’ and by contrast the remaining 11 as the 
‘poor sample’. We note that this evolution of 61200 im¬ 
plies imprints due to environment on the rich sample will 
disappear at early epochs (2 > 1 ), before they were out¬ 
liers. 

In the left hand column of Fig. [3] we compare the dis¬ 
tribution of three fundamental properties of the galaxies 
at z = 0: the gas mass (M gas ) and the star forming 
gas mass (Msf) within R 2 oo, an d the star formation rate 
(SFR) for material within 30kpc averaged over the last 
500 Myr. These quantities are all plotted as a function of 
the present-day stellar masses (M star ) and we addition¬ 
ally mark the galaxies with stellar discs. 

The cosmological trend of M gas to rise with M star is 
visible even within this restricted stellar mass sample, 
and indeed the significance of the trend shown is 4.4<r. 
For Mgp and SFR the trends are insignificant, partly due 
to the higher stochasticity of these quantities (note the 
dynamic range of Msf and the SFR is much higher than 
just for the gas) that obscures the trend in such a limited 
sample, both in terms of small number and restricted 
stellar mass range. More interestingly, perhaps, is that 
two of the rich sample (G1 and G2) are the outliers in 
Mg a s, Msf and SFR, and if we only consider SFR then 
the entire rich sample is extremal. This suggests the 
possibility that the richer environments of these galaxies 
are contributing to higher star formation rates than those 
found in the poor sample. There does not appear to be 
corresponding trends with the morphology (i.e. whether 
the stella r com ponent lies in a disc) as was also noted 
by |Few et al. (2012), and indeed in Scannapieco et al. 
(submitted) we find morphology to be largely determined 
by the merger history. 

These hints of an effect of the denser environment are 
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Fig. 3.— Comparison between the properties of the galaxy sam¬ 
ples at z = 0. Left column from top to bottom shows the gas mass 
(Mgas), the star forming gas mass (Msf) and the star formation 
rate (SFR), all as a function of stellar mass (M s t a r)- Right column 
shows the SFR and the specific star formation rate (SSFR) with 
^1200 respectively. Filled symbols indicate the galaxies that have 
stellar disc, diamonds are those in the rich sample and grey-shaded 
regions , where present, indicate the significant linear regressions 
with drier errors. 


borne out by classification with environmental overden- 
sity (Fig. |3]RHS). The SFR vs. overdensity has signifi¬ 
cant correlation, rejecting the null hypothesis at 3.8er. To 
account for comorbidity of 61200 with M star , we addition¬ 
ally considered the specific star formation rate (SSFR, 
for which we simply use SFR/M star ) and M star vs. 61200 
(not shown). The significance of the SSFR trend is lower 
( 2 . 6 cr), though both are more significant than the cor¬ 
relation of 6 i 2 oo-M star ( 0 . 8 er), i.e. environment is a bet¬ 
ter indicator of SFR than stellar mass. The correlation 
between M 2 oo and SFR is also poor as the virial mass 
closely traces M star . 

Having seen an effect of environment on SFR for MW- 
rnass galaxies at 2 = 0 it is interesting to consider 
whether this exists for their progenitors at higher red- 
shifts. Our simulations also provide evidence here, al¬ 
though we must be somewhat cautious as galaxies are 
correlated with their progenitors. In Fig. [4] we plot the 
SSFR vs 61200 for our 14 galaxies at 2 = 2, 1 and 0.5. 
At higher redshifts, the galaxies live in sparser environ¬ 
ments with higher SSFRs, and noticeably the 61200 -SSFR 
correlation fades towards earlier epochs. 

Comparing Figs. HE we can see that from 2 < 1 the 
rich sample transitions to a denser environment with a 
modest drop in SSFR whilst the poor sample has a more 
pronounced fall in SSFR but less evolution in 61200 - This 
suggests a possible explanation for the trends in Fig. [3] 
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Fig. 4.— As for the lower right panel of Fig.[3]but for redshifts 
2 = 2, 1, 0.5 from top to bottom. 


4. CONCLUSIONS 

In this letter we performed a hydrodynamical con¬ 
strained simulation of the structures within 10 Mpc of 
the MW to search for Mpc-scale environmental effects 
on galaxy properties. We compare the star formation 
rates, morphologies and gas fractions of 6 galaxies of 
MW-like stellar masses to 8 from unconstrained simula¬ 
tions and find systematically higher star formation rates 
in the galaxies that live in richer environments, though 
this trend did not extend to morphologies. 

Most notably the simulated MW and M31 candidates 
are overabundant both in present-day SSFR and in en¬ 
vironmental density, marking them as outliers compared 
to the remaining galaxies. These exceptional values dis¬ 
appear at earlier cosmological times, corresponding to 
their environmental assembly. The third galaxy which 
displayed this trend inhabited the dense environment of 
a filamentary structure. 

Although the number of galaxies is small and larger 
samples are needed to confirm these trends, our results 
suggest that galaxies of a given mass that live in richer 
environments could more easily replenish their gas reser¬ 
voirs enabling higher star formation rates. This is con¬ 
sistent with a picture where the morphologies of MW- 
like galaxies are primarily set by their merger history, 
yet environment still plays a role in their star formation 
histories. This demonstrates the need to understand in 
detail the effects of environment on galaxy properties if 
we are to understand the formation of our own galaxy. 


the cold gas is driven by the richness of the environment, 
i.e. the accretion of dense structures allows the star form¬ 
ing gas to replenish and higher rates of star formation to 
continue to z = 0. In our simulations subh alos deposit 


little of this due to their low gas content (Nuza et al. 


2014b), although they may trigger star-formation indi¬ 


rectly. Intriguingly, our trends for SFRs and star form¬ 
ing gas masses are in opposition to that of clusters vs. 


McGee et al. 

2008 

though see also 

_ oiz'- sr — v 

Ziparo et al. 

2013) 


pared to our sample, and in the case of the LG galaxies 
from the fact that the LG has a late assembly, and so 
environmental effects will be less pronounced compared 
to more evolved systems. 
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